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Abstract. We present an analytical model that characterizes two-photon transitions in the 
presence of autoionising states. We applied this model to interpret resonant RABITT spectra, 
and show that, as a harmonic traverses a resonance, the phase of the sideband beating 
significantly varies with photon energy. This phase variation is generally very different from 
the TT jump observed in previous works, in which the direct path contribution was negligible. 
We illustrate the possible phase profiles arising in resonant two-photon transitions with an 
intuitive geometrical representation. 


1. Introduction 

Autoionising states (AI) are hallmarks of electronic correlation, which shapes the reactivity of 
all many-body systems. AI states have been the subject of extensive investigation since Madden 
and Codling reported the asymmetric profile of helium doubly excited states in the first energy- 
resolved x-ray photoabsorption spectrum recorded using synchrotron radiation, a pioneering 
experiment which signed the birth of modern photoelectron spectroscopy [1]. Several years 
before, Fano had developed a model [2] in which he explained the asymmetric profiles in atomic 
photoelectron spectra, such as those seen by Madden and Codling, as interferences between two 
one-photon paths, a direct one, from the ground state to the continuum, and an indirect one, 
from the ground to the metastable state to the continnum. Synchrotron radiation gave access to 
the study of one-electron processes with unprecedented detail; the width and energy of several 
autoionising states have been accurately measured. Due to its typical properties (incoherent 


light pulses, with a duration of several picoseconds, that are highly monochromatized before 
impinging on the sample), however, synchrotron radiation did not provide access to the energy 
dependence of resonant transition phases. Hence, the dynamical information of the resonant 
process, which is encoded in such phases, was lost. 

The advent of attosecond light sources opened the way to monitor and control the electron 
motion in atoms and molecules at their intrinsic timescale [3]. In particular, such coherent ultra¬ 
fast sources of light provide a new means to investigate the role of AI states in atomic transitions 
since, hrst, the coherence of the pulses used permits to experimentally access the phases of the 
transitions and, second, the attosecond resolution permits to follow the temporal evolution of 
the AI states on a time scale smaller than their lifetime. The technique of reconstruction of 
attosecond beating in two-photon transitions (RABITT ) [4], which has already been successfully 
employed to study resonant transitions [5,6], is particularly indicated to this task. So far, 
however, only transitions through either electronic bound states [5] or autoionising vibronic 
states without any appreciable contribution from the intermediate continuum (direct path) [6] 
have been considered. In those conditions, the phase of the complex transition amplitude is 
expected to undergo a change of vr as a function of the detuning of the pump harmonic from 
the intermediate resonant state, and the measurements were indeed found to be compatible 
with this expectation. There remained a need to explore the case in which both the bound 
and the continuum intermediate states contributed to the two-photon transition. To do so, it 
was necessary to extend Fano’s model to the multi-photon hnite-pulse formulation required by 
modern attosecond interferometric techniques. A two-photon hnite-pulse resonant model, in 
particular, serves as a framework to interpret the phases of resonant transition amplitudes in 
current experiments. 

In this paper, we provide a derivation of the latter theoretical model, which was originally 
presented in [7], and which is able to reproduce to a great accuracy time-resolved resonant two- 
photon transitions. Indeed, the model has already been successfully used to interpret two recent 
RABITT experiments on resonant transitions in Helium [8] and Argon [9]. We also illustrate a 
geometrical construction that permits to interpret the phase prohles obtained from experiments. 
For a more complete set of results and a more detailed derivation, the reader is referred to 
Refs. [7,9,10]. 

2. Theory 

The lifetimes of autoionising states are often comparable to the duration of the APT and of the IR 
pulses used in RABITT spectroscopy. Furthermore, the central frequency of the harmonics may 
not coincide with a nominal multiple of Wir. For these reasons, to make quantitative predictions 
it is necessary to use a hnite-pulse formulation of the the perturbative transition amplitudes. In 
the present section, we provide a derivation for the time-dependent two-photon resonant model, 
in the simplihed case of only one isolated intermediate resonance. The derivation accounts also 
for a direct dipolar coupling between the localised component of the intermediate metastable 
state and the hnal continumm. 

The time-dependent lowest-order perturbative amplitude for the two-photon transition from 
an initial atomic bound state \g) to a hnal non-resonant continuum state \"fE), for a linearly 
polarized held F{t) = eF{t), is 

•^'yE,g = y dujF{F - Fg - Uj)F{u})M^E,g{‘^), (1) 

where M.yE,gi^) = {iF\OGq {iOg +(jj)0\g) is the two-photon dipole matrix element, the Fourier 
transform (FT) is dehned as F{oj) = (27r)“^/^ / F{t) exp{iujt)dt, Gq{F) = {F — F[ + i0+)“^ 
is the retarded resolvent of the held-free hamiltonian Ff, O = e ■ O is the dipole operator (in 


velocity gauge, F is the vector potential and O = aP, where a is the fine-structure constant). 
The index 7 indentifies collectively all the quantum numbers, other than energy, that uniquely 
identify the final continuum, and which are needed to differentiate it from other continuum 
states, such as those populated by one-photon transition from the ground state. Atomic units 
are used throughout, unless otherwise stated. We assume that the field of the impinging 
pulses can be expressed as a linear combination of Gaussian pulses, F{t) = Foexp[—cj^(f — 
fo)^/ 2 ] cos[u;o(t — to) + </>]) where Fq, wq, to, a and (j) are the amplitude, carrier frequency, center, 
spectral width and carrier-envelope phase of the pulse, respectively. The FT of such pulse is 
F{uj) = F~^ (uj) + F~ (co), where F^{u}) = Fo(2cj ) ^ exp(iu;to) exp[—(a;=Fa;o)^/( 2 cj^)] are 

the components responsible for photon absorption and emission, respectively. In a pump-probe 
experiment with an XUV attosecond-pulse train (APT) in association with an isolated IR pulse, 
the APT center txuv conventionally defines the time origin, txuv = 0, while the center of the 
IR pulse coincides with the pump-probe time delay r, fiR = r. In the present context, the total 
external field is conveniently expressed in terms of synchronized pulses for the odd harmonics, 
which give rise to the APT, and of a delayed pulse for the IR field, F{t, r) = FaptC^) + Fi^{t—T), 
where FAPT(t) = Fii 2 n+iit)- In the RABITT scheme, where one XUV photon is absorbed and 
one IR photon is either absorbed or emitted, the transition amplitude corresponding to a given 
sideband SB 2 n is given by the sum of four contributions, each associated to a time-ordered 
perturbative diagram. Since the contribution of the two diagrams in which the IR photon is 
exchanged first is generally small and can be neglected, so the transition amplitudes become 
= —i f du;Fj^(E — Eg — uj)F^^^^^{Ljj)M.^E,g{'-jj) ■ The energy-resolved intensity of the 
sideband, Isb, is computed from the transition amplitudes as 
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To evaluate the two-photon ionization matrix element Ai^E,g(<^), we assume that the 
intermediate eigenstates of H are well represented by a single resonant channel, Hl'ijjaE) = 
IV’oe)F that can be expressed, following Fano’s formalism [2], in terms of known bound state, 

I a), and featureless continuum states, \ae), which are eigenstates of a reference hamiltonian Hq, 
Ho\a) = Ea\a), Ho\aE) = \aE)e, {aE\aE') = 6{E — E'), 
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where 14,= W\H — Ho\aE), Ea = Ea — iTa/2 is the complex resonance energy, Ea = 
Ea + P J d£\Va,ae\^/{E — e), F^ = 27r\Va^aE\‘^ [2]. For simplicity, in the present derivation we 
disregard the role played by virtual excitations of intermediate bound states. The transition 
matrix elements between a localised state and Fano continuum can be parametrized as 
{'ijjaE\0\g) = OaE,g{eE + q){eE-i), where eg = 2{E-Ea)/Ta and q = Oag/{EVaEOEg) [2]. The 
two-photon ionization matrix element, therefore, can be written as 
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By applying the on shell approximation, {'^E\0\a£) ~ 0^a{E)5{E—£), which, for singly-charged 
parent ions, is quite accurate at energies of the order of 1 a.u. above the threshold or larger, 
and assuming that both O^aiE) and Va^aE are sufficiently slowly varying functions of E, it is 
immediate to see that 
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We can now insert this expression in the two-photon matrix element, 
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The integral in this last expression can be easily computed closing the integration circuit in the 
lower half of the complex plane and applying Cauchy’s residual theorem, 
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ag 


where ojij = oji — uij and we introduced the parameter (3a = T^O^E,aya,oLEl0^a. 

If one is interested in the long-pulse limit only, for which the harmonic spectrum is strongly 
peaked at Wh, the resonant and non-resonant two-photon transition matrix elements, for the 
absorption of an XUV photon Wh followed by the absorption/emission of one IR photon, cuia, 
can be approximated with their value at w = Wh > 
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where we used the energy-conservation principle E = Eg E Wja and we introduced the 
notation Ch = eEg+ui^i reduced detuning of the harmonic from the resonance. Notice 

that when the intermediate autoionizing state |o) is not directly radiatively coupled to the final 
continuum, i.e., /3a = 0, we recover an expression equivalent to Eq. (4) in the main manuscript, 


A4^ (oj) ~ (nE q 


( 9 ) 


The dimensionless parameter (a = 2a;ia/3a/ra = O^E,a/((^'ya(^m^yaa) expresses the relative 
strength of the direct |a) —)• \"fE) radiative transition compared to the indirect path for the same 
process |a) —)• \aE) —>• lyE), in which the system first decays non-radiatively to the continuum 
and subsequently exchanges a photon in a continuum-continuum transition. As Ch increases 
from —oo to -|-oo (i.e., from large negative to large positive detuning of the resonant harmonic) 
the resonant factor (ch -|- g^)/(eH E i) describes, counterclockwise, a circular trajectory in the 
complex plane, which starts from (1,0), is centered at (1 — zg^)/2 and has radius |1 — iq^\/2, 
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The circle intercepts the origin only if ^a = 0. If Ca ^ 0, the origin falls outside the circle and 
hence the phase of the resonant matrix element g experiences a continuous excursion with 
no net variation, while if Ca^O, the circle encloses the origin, and the phase of g undergoes 
a smooth overall excursion of 27r. In the long-pulse limit, therefore, when ((a 7^ 0, if the phase 
of the resonant amplitude for the absorption of one IR photon performs a jump of 27r, that for 
the emission of an IR photon will experience no net variation, and viceversa. 

With finite pulses, the transition amplitude is given by the convolution (1) of M.^E,g{(o) with 
the FTs of the held. The calculation for Gaussian pulses is lengthy but straightforward, and 
the result can be expressed in closed form in terms of the Faddeeva special function, w{z), here 
dehned as the analytic continuation of ztt”^ f^^dtexp{—t‘^)/{z — t) {^z >0). In the case of 
the absorption of an XUV photon Wh followed by the absorption (-I-) or emission (—) of an IR 
photon WiR, the result is 
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Figure 1. Trajectory of the resonant two-photon amplitude in the complex plane (upper panels) 
and corresponding phase (lower panels), as the reduced detuning e = 2{E — Ea)/Ta increases 
from large negative to large positive values. In the left panels, the only contribution to the 
transition comes from the intermediate bound component {q = oo). In the central panels, both 
the intermediate continuum and bound components contribute to the dipolar transition from the 
ground state (g = 1), but the intermediate bound component is not radiatively coupled to the 
final continuum = 0. Finally, in the right panels, this final restriction is removed = 0.2), 
and, as a consequence, the circle of the transition amplitude to the upper sideband is contracted 
and misses the origin, thus experiencing no net phase change, while the circle of the transition 
amplitude to the lower sideband gets expanded, thus encircling the origin and experiencing a 
net phase transition as large as 27r. See text for more details. 


where J-{t) is an inessential form factor of the field, and 


± ^ _i_ M [E - (Eg+ oj^±ojjji)]/al + iT / _2 , I 2 
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It is instructive to verify that, in the limit of long {atTa S> 1) overlapping (r <C at) pulses, 
we recover expressions similar to the ones seen above, as expected. Indeed, in this limit, 
the argument of the Faddeeva function diverges, so one can use the asymptotic expansion 
w{z) ~ . At the nominal centre of the two-photon signal {E = Eg + ui-n E Wm), 

therefore, zf^ ~ at {Eg + Uyi — e) /\/2. Using the energy-preserving condition, one can simplify 
the transition amplitudes to 


(13) 

With finite pulses, the resonant amplitudes describe contracted circles. For sufficiently 

small values of |Ca|) therefore, both amplitudes would fail to enclose the origin, and hence the 
jump of 27r wouldn’t be observable. 

3. Conclusions 

We have derived a model which permits us to predict the phase variation of two-photon resonant 
amplitudes, extracted from current attosecond pump-probe experiments. In the limit of long 
pulses, we find an expression similar in form to that for the resonant one-photon transitions, but 
which requires an additional parameter C, arising from the exchange of an IR photon with the 
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Figure 2. Sideband photoelectron signals in a RABITT pump-probe photoionization of the 
helium atom from the ground state, as a function of both the photoelectron energy (vertical 
axis) and the pump-probe time delay (horizontal axis). Each panel corresponds to a different 
reduced detuning Ch of harmonic H41 from the 2s2p intermediate doubly excited state 
{E = 35.56 eV, E = 0.037 eV, q = —2.77, C/'^m = 0-19 eV“^), whose energy is indicated by 
the thin white dashed line: (a) Ch = —43, (b) Ch = —25, (c) Ch = 0.9, (d) Ch = 17, (e) Ch = 35. 
At large negative (a) or positive (f) detunings, the upper and lower sidebands are in phase to 
a very good approximation. As the harmonic traverses the resonance, on the other hand, the 
beating of the upper and lower sidebands get clearly out of phase and the sideband profile itself 
is distorted. These panels are computed with an extension of the model described in the text, 
which accounts for final resonant states. The thin horizontal feature above 37 eV, visible in 
panels (c-e) is due to the final 2p^ autoionising state. 


autoionising state. The model has a compact analytical finite-pulse formulation, which has been 
instrumental to interpret recent findings in resonant RABITT experiments [8,9]. With the use 
of the model and simple geometrical constructions, we have characterised the possible profiles 
of the amplitude phase that can be expected in resonant two-photon transitions as a function 
of detuning. These can vary from a finite excursion with no net phase variation to a total jump 
of 27r, depending on the relative strength and phases of the atomic transitions involved, the 
duration of the pulses used and the number of intermediate open channels. 
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